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In this paper, we have considered a model of modified Chaplygin gas in VSL theory with variable 
gravitational constant G. We have shown that the evolution of the universe starts from radiation 
era to phantom model. The whole evolution of the universe has been shown diagramatically by 
using statefinder parameters. 

^ ' PACS numbers: 98.80.Cq, 95.35.+d 
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(N ! I. INTRODUCTION 

a : 

Recently, the matter field can give rise to an accelerated expansion [1] for the universe stems from the 
observational data regarding the luminosity-redshift relation of type la supernovae [2,3] up to about z ~ 1. This 
matter field is called quintessence matter (shortly Q-matter). This Q-matter can behave like a cosmological 
constant [4,5] by combining positive energy density and negative pressure. So there must be this Q-matter 
either neglected or unknown responsible for this accelerated universe. At the present epoch, a lot of works 
■ has been done to solve this quintessence problem and most popular candidates for Q-matter has so far been 
' a scalar field having a potential which generates a sufficient negative pressure. The Q-matter behaves as a 
perfect fluid with a barotropic equation of state and so some effort has been investigated in determining its 
adiabatic index at the present epoch [6, 7]. Chimento et al [8] have showen that a combination of dissipativc 
effects such as a bulk viscous stress and a quintessence scalar field gives an accelerated expansion for an open 
universe (k = —1). Banerjee et al [9] also have shown that it is possible to have an accelerated universe with 
BD-theory in Friedmann model without any matter. 
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Another alternative candidate for Q-matter is exotic type of fluid — the so-called Chaplygin gas which obeys 
5h ■ the equation of state p = —B/p, (B > 0) [10], where p and p are respectively the pressure and energy density. 
Subsequently, the above equation was generalized to the form p — —B/p a ,0 < a < 1 [11, 12] and recently it 
, was modified to the form p = Ap — B/p a ,(A > 0) [13, 14], which is known as Modified Chaplygin Gas. This 
model represents the evolution of the universe starting from the radiation era to the KCDM model. 

The possibility that the speed of light c might vary has recently attracted considerable attention. In a 
cosmological setting, the variations in c have been shown to solve the cosmological puzzles - the horizon, 
flatness and Lambda problems of big-bang cosmology. The variations of velocity of light can also solve the 
quasi-lambda problems. For power-law variations in the velocity of light with the cosmological scale factors, 
Barrow et al [15] have shown that flatness problem can be solved. The Machian VSL scenario in which 
c = cod™, introduced by Barrow [16] has significant advantages to the phase transition scenario in which the 
speed of light changes suddenly from c~ to c + , preferred by Albrecht and Magueijo [17]. For changing c, the 
geometry of the universe is not affected. In classical electromagnetism the speed of light is only constant 
in vacuum and it 'varies' in dielectric media. The changing c breaks Lorentz invariance and conservation of 
energy. The geometry of Universe is not affected by a changing c. One can allowed a changing c to do the job 
normally done by superluminal expansion. For changing c, the gravitational laws should be modified. The basic 
assumption is that a variable c does not induce corrections to curvature in the cosmological frame and therefore, 
Einstein's equations, relating curvature to stress energy are still valid. The reason behind this postulate is 
that c changes in the Local Lorentzian frames associated with cosmological expansion. The effect is a special 
rclativistic effect and not a gravitational effect. Therefore curvature should not be related with the variation of c. 

The organization of the paper is as follows: In section II, we have studied field equations and its solutions 
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due to VSL theory for modified Chaplygin gas with varying G. In section III, we have discussed whole stage of 
the evolution of the universe using statefmder parameters. In the last section (i.e., section IV) we give some 
remarks of this paper. 



II. BASIC EQUATIONS AND SOLUTIONS 

In VSL theory and varying G, the Einstein equations for FRW model is 



a2 kc 2 (t) _ 8irG(t)p 



and 

a _ AnG(t) ( 3p 



c 2 (t) 



(2) 



where p and p are the pressure and energy density respectively. Here we consider modified Chaplygin gas in 
VSL theory as 

p = Apc 2 (t) - {p f^ }a , A > 0, B > 0, < a < 1. (3) 
The energy conservation equation incorporating time variation in c(t) and G(t) is 

„a / P\ G 3kcc 

Now assume, the velocity of light and G are in power-law form of scale factor as 

c(t) = c a n (5) 

and 

G = G a m (6) 

where cq and Go are both positive constants. Since we know that speed of light decreases with time and G 
increases with time, so that n must be negative and m must be positive. 

Using equations (3), (4), (5) and (6) we have the solution of p as (for k = 0) 



P = 



(7) 



where M = — , n , n) , N = -Sj and ci is a constant. 

(3+3A-2n+m)c 2 {1+a> ' S„ 

Putting the value of p from (7) in (2), we have 

a 2 da 

c 2 t (8) 
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[M a ( 3 - 2n + m + 3A )( 1 +«) + Afj^+sy 

where C2 = \J~- 



'8ttG 
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The solution of scale factor a(t) has the form 
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„ffi±^ F 3(1 + A) i , 3(1 + A) M n 3+3A+ m -2 n] _ 3 n , , V NWT^J t 

(9) 

where 2 -F\ is the hypergeometric function. 

Now for small value of scale factor a(t) i.e., for early universe, p is very large and corresponds to the universe 
dominated by an equation of state p = Ape 2 . 

Also for large value of scale factor a(t) i.e., for late universe, the relation between p and p satisfies the equation 
of state 

P=-\l + -^—) P c ( 10 ) 

So evolution of the universe starts from radiation era (A = 1/3) to phantom model. From observations, we 
can conclude that m < 2(n + 1). For ACDM model we have m = 2n. 



III. STATEFINDER PARAMETERS AND EVOLUTION OF THE UNIVERSE 

In 2003, V. Sahni ct al [18] proposed a pair of parameters {r, s}, called statefinder parameters. In fact, 
trajectories in the {r, s} plane corresponding to different cosmological models demonstrate qualitatively different 
behaviour. The above statefinder diagnostic pair has the following form: 

where H (= -) and q (= — £§) are the Hubble parameter and the deceleration parameter respectively. These 
parameters allow us to characterize the properties of dark energy. Trajectories in the {r, s} plane corresponding 
to different cosmological models, for example ACDM model diagrams correspond to the fixed point s = 0, r = 1. 

The expressions for q, r and s are in the following: 

1 + 3A 3B 



2(pc 2 ) 1 +« 

1 + 3A 3B 



(12) 
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(14) 



The expressions of r and s are very complicated, so we can not found the relation between r and s in closed 
form. 
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Fig. 1 shows the variation of q against a for A = 1/3, a = 1/2, B = 1, n = —1/2, m = 1. The plot of g 
against a reveals that the deceleration parameter indeed has a sign flip in the desired direction and indicates 
an early deceleration (q > 0) followed by a late time acceleration (q < 0) of the universe. From fig. 2 we see 
that the curve starts from the radiation era and goes asymptotically to he dust model for s > and r > 7. 
But for r < 7, the portion of the curve represents the evolution from dust state (s = — oo) to phantom model 
(s > 0) and this curve crosses ACDM stage (r = 1, s = 0). Thus the total curve represents the evolution of the 
universe starting from radiation era to phantom model. 
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IV. DISCUSSIONS 

In this work, we have proposed a model for modified Chaplygin gas in VSL theory with variable G. 
In this model, we have found the evolution of the universe starts from radiation era to phantom model. 
If m = n = i.e., c =constant and G ^constant then we can recover our previous results [14]. Thus 
variable c and G can drive the evolution of the universe from radiation to phantom stages instead of ACDM 
stage. From equation (10) we see that for ACDM model m — 2n and for phantom model m < 2(n + 1). 
If G ^constant or c =constant then any one of variable G and c can drive the evolution of the universe 
into the phantom model. Thus variable G or c has significant role for the evolution of the universe in large stages. 
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